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One-way surface magnetoplasmon cavity and its application for nonreciprocal devices 
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We theoretically analyze surface magnetoplasmon modes in a compact circular cavity made of 
magneto-optical material under a static magnetic field. Such a cavity provides two different physical 
mechanisms for the surface wave to circulate in a unidirectional manner around the cavity, which 
offers more freedom to realize one-way surface wave. We also show the interaction between this 
one-way cavity and waveguides, through an example of a circulator, which lays the fundamental 
groundwork for potential nonreciprocal devices. 
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Surface plasmons (SPs) at metal-dielectric interfaces 
enable confining and manipulating electromagnetic waves 
on a subwavelength scale [l|. A variety of applications 
based on the concept of SPs have been achieved at dif¬ 
ferent frequency regimes such as focusing of light beyond 
the diffraction limit @ , ultracompact nano-lasers [3] and 
nano-antennas 0 , 0 . Applying a static magnetic field 
on SPs breaks the time-reversal symmetry and gives rise 
to the nonreciprocal propagation of SPs called surface 
magnetoplasmons (SMPs) pi8[. The asymptotic fre¬ 
quencies of SMPs in the forward and backward direc¬ 
tions are split by the external magnetic field. This pro¬ 
vides a one-way regime for SMPs, i.e., the propagation of 
the surface wave is only allowed in one direction within 
the frequency range between two asymptotic frequencies 
UM- There are also other mechanisms supporting one¬ 
way waves [IM3 including topological protected edge 
states in photonic crystal. Compared with these however, 
the mechanism based on SMPs seems more attractive due 
to its simple configuration. 

Recent research has shown that SMPs are immune 
to backscattering from disorder 0 Ezl. which might 
be useful for applications involving isolators. Most re¬ 
sults are based on SMPs at a planar interface and use 
only a single frequency range 0 G3 [IM]. Since the 
cavity structure is always compact and at the heart 
of many photonic components (e.g. High-Q surface- 
plasmon-polariton whispering-gallery microcavity 0 ), 
it is interesting to study the behavior of SMPs in cavity 
structures and investigate the physical mechanisms for 
a one-way SMP cavity mode, i.e., when SMPs circulate 
around the cavity in a unidiretional manner. Integrating 
such SMP structures with conventional waveguides is also 
of particular interest for constructing some nonreciprocal 
devices. 

In this letter, we theoretically analyze SMP modes 
in a compact circular cavity made of magneto-optical 
(MO) material under a static magnetic field and demon¬ 
strate two different physical mechanisms for one-way cav¬ 


ity modes. The first one results from splitting the asymp¬ 
totic frequencies for clockwise and anti-clockwise cavity 
modes, which forms the upper one-way frequency range. 
The other is obtained by different cut-off frequencies for 
the lowest clockwise and anti-clockwise cavity modes, 
which forms the lower one-way frequency range. The 
circulation directions of SMPs in the cavity are oppo¬ 
site in these two different mechanisms. We also study 
the interaction between such a one-way SMP cavity and 
conventional waveguides, which is fundamental for de¬ 
signing nonreciprocal devices. A compact circulator is 
designed and shown by direct numerical simulations as 
an example. 

We firstly analyze the cavity modes for SMPs to get 
an insight into the possible one-way regimes. The cav¬ 
ity structure is a circle with radius R surrounded by air 
in a two-dimensional (2D) system as shown in the insert 
of Fig. 1(a). The circular cavity is made of a semicon¬ 
ductor which is an MO material at THz frequencies. By 
applying a static magnetic field in the +z direction and 
assuming the material is lossless, the relative permittiv¬ 
ity tensor of the semiconductor takes the following form 
in cylindrical coordinates (p,p,z) [6]: 
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s oo (1 — ^1), where co is the angular frequency, uj v is the 
plasma frequency of the semiconductor, uj c = eB/m* is 
the electron cyclotron frequency (e and m* are, respec¬ 
tively, the charge and effective mass of the electron, and 
B is the applied magnetic field), and £oo is the high- 
frequency permittivity of the semiconductor. We only 
consider the TM mode (H z , E p , E p ) in this cavity, where 
E z = H p = Hp = 0. Maxwell’s equations yield 
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Taking the field H z = 'ip(p)e irnip , where m is an integer 
indicating the azimuthal mode number, we obtain the 
following differential equation from Eq. (2): 


[^2+-^- + (k 2 -^M(p) = 0 (3) 
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where k 2 = oj 2 po£o£ v , with Voigt permittivity £ v = £\ — 
£ 2 /£\. Inside the circle (p < R), the solution for Eq. (3) 
is the Bessel function of the first kind J m (kp). Outside 
the circle (p > R ), for outgoing waves, the solution for 
Eq. (3) is the Hankel function of the first kind Hm\kop ), 
where ko is the wave vector in air. Therefore, the solution 
is given by 




A * Jm(kp) P ^ R 
C ■ Hm\kop) p > R 


( 4 ) 


According to the boundary conditions that H z and E p 
are continuous at p = R, we obtain the following eigen- 
frequency equation of the cavity modes: 

JmjkR ) _ (£ 2 m/R) ■ JrnjkR ) + £1 k • J / m (fc-R) 

Hm\k 0 R) (£? - 4)ko • H$'(koR) 

We can solve Eq. (5) numerically for each given az¬ 
imuthal mode number m and obtain the dispersion re¬ 
lation between the eigenfrequency and the mode number 
m. The linear term (£ 2 m/R) • J m (kR) with respect to 
m, which originates from the off-diagonal element of the 
permittivity tensor, breaks the left-right symmetry of the 
dispersion relation. Therefore, the solutions for Eq. (5) 
and also the field distributions are different for the cavity 
modes | +m) and | —ra), where |+ra) and | —m) represent 
the modes with positive and negative mode numbers, re¬ 
spectively. In addition, we only investigate the funda¬ 
mental radial mode, since it is the most confined of the 
radial modes. 

We take a semiconductor disk made of InSb with R = 
250 p m as an example to show the solution for Eq. (5). 
At room temperature, the semiconductor InSb takes the 
following parameters: £qo = 15.6, f p = uj p /2tt = 2 THz 
and ra* = 0.014rao (mo is the free electron mass) (HJ. 
We study three cases with increasing applied magnetic 
field B = 0 T (o+ = 0), B = 0.025 T (cj c = cap/40) and 
B — 0.1 T (cj c = cjp/10) along the -hz axis and solve Eq. 
(5) numerically. The solutions for those three cases are 
illustrated in Fig. 1(a), where the dots show the disper¬ 
sion relations and the gray region represents the region 
above the air light line. The air light line is determined by 
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FIG. 1: (a) The dispersion relation between the eigenfre¬ 

quency and azimuthal mode number m of the circular cav¬ 
ity. The circular cavity with radius R is surrounded by air in 
the two-dimensional (2D) system as shown in the insert. The 
cavity is made of the semiconductor InSb, which is a magneto¬ 
optical (MO) material at THz frequencies. Three cases with 
increasing electron cyclotron frequency uj c = 0 (black square 
dots), uj c — cjp/40 (red triangular dots) and u c = cj p /10 (blue 
circular dots) are plotted. The gray region represents the re¬ 
gion above the light line. The distributions of the modal field 
H z are shown in (b) m — +17, (c) m — +30, (d) m — —17 
and (e) m — —30, respectively, when uj c — oj p /1Q. 


f = mc/2irR. The dispersion relations start from differ¬ 
ent cutoff frequencies and approach different asymptotic 
frequencies in these three cases. Firstly, in the absence 
of an external magnetic field (B = 0 T, uj c = 0, black 
square dots), the dispersion relation is symmetric, as the 
time-reversal symmetry is preserved. When m -+ ± 00 , 
the mode has an infinite azimuthal wave vector, and the 
asymptotic frequency is f s = (u p /2n)\fe^£JJho£+££ir ), 
which is the same result for SPs without magnetization. 
When m is very small, due to the limitation of the air 
light line, the corresponding eigenfrequency can be a 
complex number, which corresponds to a radiative mode 
:23]. The lowest mode with a real eigenfrequency gives 
the cutoff frequency for nonradiative modes. In this let¬ 
ter, we only consider the nonradiative modes. Note that 
when R — y 00 , the cutoff frequency will approach zero, 
which is exactly the case for SPs at a planar interface. 
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The curvature of the interface enables one to achieve the 
cutoff frequency above zero. For the second case with 
B = 0.025T (cj c =8 cjp/40, red triangular dots), the dis¬ 
persion relation is asymmetric, and the asymptotic fre¬ 
quencies are split by the external magnetic field. The 
+m branch of the dispersion relation is higher than the 
dispersion relation with B — 0 T, while the — m branch 
is lower than the one with B = 0 T. When m —>• Too, ap¬ 
plying the formulas J m _^ +00 (,s) ~ (l/ \[2jrm) • (es/2ra) m , 

h£\ +00 (s) ~ {-i^/2/y/mn) ■ ( es/2m )~ m , J- m (s) = 
(-1 ) m J m {s) and H- m (s) = (-1 ) m H m (s) in Eq. (5), 
we obtain the following asymptotic frequencies 


which are the same as SMPs at a planar interface dis¬ 
cussed in 0. Within the one-way frequency range 
/_oo < / < /+ oo, there only exists \+m) modes, and 
SMPs can only circulate around the cavity anti-clockwise. 
Furthermore, there does exist different cutoff frequencies 
for |+ra) and | —m) cavity modes. The cutoff frequency 
f+ c of the |+ra) mode is higher than the cutoff frequency 
f- c of the | —m) mode, so there is another one-way fre¬ 
quency range. Within /_ c < / < /+ c , there only exists 
| —m) modes, and SMPs can only circulate around the 
cavity clockwise. In the third case, we apply a stronger 
magnetic field B = 0.1 T (cj c = cj p / 10, blue circular 
dots). The applied strong magnetic field considerably 
breaks the symmetry of the dispersion relation. The cut¬ 
off frequency f+ c of the |+m) mode can be lifted higher 
than the asymptotic frequency /_^ of the | —m) mode. 
In this case, the modes \+m) and \—m) are separated in 
completely different frequency ranges. The |+m) modes 
are within / +c < / < /+ oo and SMPs can only circulate 
anti-clockwise in this frequency range. The \—m) modes 
are within /_ c < / < /_^ and SMPs can only circu¬ 
late clockwise in this frequency range. The circulation 
directions are opposite in these two different one-way fre¬ 
quency ranges. This interesting feature offers more free¬ 
dom, both in the one-way direction and in the frequency 
range, for designing nonreciprocal components. 

To verify the theoretical results obtained from Eq. (5), 
we also solve the eigenfrequency of the cavity mode and 
calculate the modal field distributions by a finite element 
method (FEM) in the commercial software COMSOL. 
The deviation of the FEM results from the results ob¬ 
tained by Eq. (5) is less than l%o- The modal field 
(H z ) distributions of the cavity modes |d=17) and |±30) 
for B = 0.1 T are plotted in Fig. 1(b)-1(d), respectively. 
The field is confined at the boundary of the cavity and 
this confinement increases with \m\. We also note that 
the field distributions are different for the | -bra) and | —m) 
modes, since the strong external magnetic field breaks 
the symmetry. To clearly demonstrate the circulation di¬ 
rection of SMPs, we made videos of the harmonic wave 
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FIG. 2: The schematic of the circulator. The circulator con¬ 
sists of a circular cavity with R = 250 um and two straight 
dielectric waveguides. The straight waveguides are symmet¬ 
ric about the center of the cavity. The width of the straight 
waveguides is W = 40 um, the length is L = SR and the 
dielectric constant of the waveguides is Sd = 5. In addition, 
the gap between the cavity and the waveguide is g = 20 um. 


propagation around the cavity at the eigenfrequencies of 
modes |d=17). To excite these modes, a magnetic cur¬ 
rent point source is placed at the boundary of the cavity. 
The field H z of the wave is recorded (see supplementary 
material online). 

In this section, we study the interaction between the 
one-way SMP cavity and the conventional waveguides, 
which is fundamental for designing nonreciprocal devices. 
We use the one-way SMP cavity to construct a circula¬ 
tor as an example. The circulator consists of a circular 
cavity and two straight dielectric waveguides as shown 
in Fig. 2. The parameters of the cavity are the same 
as those in the third case discussed above. The straight 
waveguides are symmetric about the center of the cavity. 
We denote the width, length and relative permittivity of 
the waveguides by IF, L and e<i- The gap between the 
cavity and the waveguide is g. The dispersion relation 
for the fundamental guiding TM mode (H z , E x , E y ) in 
the waveguide can be obtained by solving 

fci ,k{W 

k 2 =— tan(—-) (7) 

£d A 

where fci =A/u) 2 Mn£n£d — /3 2 , fc^=\//3 2 —u; 2 Mn£n and /3 is 
the propagation constant [24| . To achieve a strong inter¬ 
action between the guiding modes and the cavity mode, 
the phase matching condition |2fj] is required: 

P(fm) = Pe(m)(l ~ (8) 

where f m is the eigenfrequency of the cavity mode |m), 
and /3 e (m) is the equivalent propagation constant of the 
cavity mode | m). When the modal field is confined 
well at the boundary of the cavity, /3 e =m/R. The term 
(l—g/2R) includes the effect of the curvature at the cav¬ 
ity boundary on the propagation constant as seen by the 











4 





TTi 

w 

O 

(a) 


(b) 




O 

A-, 

0, 

(d) 


(c) 





Frequency (THz) 


FIG. 3: Simulation result of wave propagation at the reso¬ 
nance frequency for m = +17. In (a)-(d), input power is at 
Port 1, 2, 3 and 4, respectively. 



FIG. 4: Simulation result of wave propagation at the reso¬ 
nance frequency for m = —17. In (a)-(d), input power is at 
Port 1, 2, 3 and 4, respectively. 
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FIG. 5: Calculated S parameters for the circulator, (a) S 
parameters within the upper one-way frequency range. The 
dips are the resonance frequencies corresponding to m — 
+ 16,+17,+18. (b) S parameters within the lower one-way 
frequency range. The dips are the resonance frequencies cor¬ 
responding tom = —16, —17, —18. 


straight waveguide approximately 0. Due to the asym¬ 
metric dispersion relation of the cavity mode, coupling 
occurs in a uni-directional manner between the guiding 
mode and the cavity mode within the one-way frequency 
ranges. We set m = 17 and find that the structure with 
W = 40 /im, Ed= 5 and g = 20 /im satisfies Eq. (8). In 
addition, the length is L = SR = 2000 gm and Ports 1-4 
are labeled at the end of the waveguides. 

Next, we study the performance of the circulator at 
both one-way frequency ranges. As the radius of the 
cavity is in the order of wavelength, and the structure 
of the circulator is compact, we use the FEM method 
to simulate the performance of the circulator. For the 
upper one-way frequency range, in which only |+m) ex¬ 
ists, Fig. 3(a)-3(d) show the wave propagation at the 
resonance frequency for m = +17, when the input power 
is at Port 1, 2, 3 and 4, respectively. The input wave 
from Port 1 is not coupled with the cavity mode and 
goes out at Port 2 [Fig. 3(a)], while the input wave 
from Port 2 is coupled with the cavity mode and goes 
out at Port 3 [Fig. 3(b)]. The wave travels in the or¬ 
der Portl+>Port2+>Port3—>Port4-+Portl. For the lower 
one-way frequency range, in which only \—m) exists, Fig. 
4(a)-4(d) show the wave propagation at the resonance fre¬ 
quency for m = —17, when the input power is at Portl, 
2, 3 and 4, respectively. The wave travels in the order 
Portl+>Port4+>Port3->Port2+>Portl. 


The properties of transmission and isolation are shown 
by the calculated S parameters in Fig. 5. Within the 
upper one-way frequency range [Fig. 5(a)], the transmis¬ 
sion coefficient S 21 (red solid line) from Port 1 to Port 
2 is close to 0 dB, which indicates that the guided wave 
from Port 1 is not coupled with SMPs in the cavity. The 
small losses (small dips of S 21 ) are caused by a weak 
coupling with some other whispering-gallery-like modes 
inside the cavity. This coupling can be eliminated by in¬ 
creasing the gap between the cavity and the waveguides 
(not shown here). The three dips of the isolation S 12 
(blue dashed line) are the resonance frequencies corre¬ 
sponding tom = +16, +17, +18. The resonance frequen¬ 
cies are very close to the eigenfrequencies of the cavity 
modes |+16), | + 17), |+18), as the cavity modes are per¬ 
turbed slightly by the 2 straight waveguides. Within the 
lower one-way frequency range [Fig. 5(b)], the transmis¬ 
sion coefficient 612 (blue dashed line) from Port 2 to Port 
1 is close to 0 dB. The three dips of the isolation £21 (red 
solid line) are the resonance frequencies corresponding 
to m — —16, —17, —18. The resonance frequencies are 
shifted a little from their corresponding eigenfrequencies 
of | —16), | —17) and | —18), since the geometric structure 
is designed to fulfill the phase matching condition Eq. 
( 8 ) at m = +17. All the dips are lower than -20 dB, 
indicating good isolation performance. 

In conclusion, we have theoretically analyzed the SMP 
mode in a compact circular MO cavity under a mag- 
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netic field. In such a cavity, we have found that the 
different asymptotic and different cut-off frequencies for 
clockwise and anti-clockwise modes lead to two one-way 
frequency ranges for SMPs in the cavity. These multiple 
mechanisms for achieving one-way SMPs offer more free¬ 
dom, both in the one-way direction and in the frequency 
range, for designing nonreciprocal photonic components. 
We also have studied the application of this cavity in a 
four-port circulator as an example to show the interac¬ 
tion between the one-way SMP cavity and waveguides, 
which is the heart of nonreciprocal devices. We believe 
that the mechanisms we found in the cavity will provide 
people more ways to manipulate SPs and our idea may 
be useful in a variety of potential applications ranging 
from THz signal isolation to isolators on chips. 
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